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The interaction of high energy cosmic rays with the Earth's atmosphere produces extensive air showers of secondary 
particles with a large muon component. By exploiting the sensitivity of neutrino telescopes to high energy muons, it 
is possible to use these detectors for precision cosmic ray studies. The high rate of cosmic-ray muon events provides a 
high-statistics data sample that can be used to look for anisotropy in the arrival directions of the parent particles at the 
per-mille level. This paper reports on the observation of anisotropy in the cosmic ray data collected with the IceCube 
neutrino telescope in the 20-400 TeV energy range at multiple angular scales. New data from the IceTop air shower 
array, located on the ice surface above IceCube, shows an anisotropy that is consistent with the high-energy IceCube 
results. The sensitivity of IceTop to all the components of the extensive air shower will allow us to explore in more detail 
the characteristics of the primary cosmic rays associated with the observed anisotropy. 
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It has been almost 100 years since the discovery of cos- 
mic rays by Victor Hess. The origin of these energetic 
particles, however, still remains an enduring problem in 
astrophysics. Based on indirect evidence, it is believed 
that cosmic rays (CRs) with energies up to a few PeV 
are accelerated in supernova remnants distributed accross 
our galaxy. A direct test of this hypothesis is a challenge 
since the arrival directions of cosmic rays at Earth do not 
point back to their sources due to the scrambling action 
of the galactic magnetic field (GMF) during propagation. 
For this reason, the search for the sources of cosmic rays 
is usually performed either at the highest energies where 
the influence of the GMF is small, or by making use of a 
neutral messenger particle as the neutrino. 

Even if a direct detection of cosmic ray sources is not 
feasible at or below PeV energies, their discrete spatial dis- 
tribution should create an observable dipolar anisotropy of 
per-mille strength pQ [2] [3] . The energy dependence of the 
phase and amplitude of this kind of anisotropy would be 
dominated by details in the propagation process such as 
the geometry of the galaxy, the energy dependence of the 
CR diffusion coefficient, and the age and injection spec- 
trum of the sources. Other factors, such as turbulent prop- 
agation in the GMF [4], heliospheric effects [5 , and spe- 
cial magnetic field configurations [6] [7], could give rise to 
anisotropy at smaller angular scales. A different process, 
known as the Compton-Getting effect [8], could also create 
a dipole anisotropy due to the relative motion of the solar 
system with respect to the cosmic ray plasma. 
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Several experiments in the northern hemisphere have 
reported on the observation of anisotropy at TeV energies 
[9] [10] [11] [12] [13] . Two features dominate the northern sky 
in cosmic rays: a dipole-like large scale structure with an 
amplitude of ~ 10 -3 , and a small scale anisotropy with 
significant structure at angular sizes between 10° and 30°. 
The observed dipole anisotropy is inconsistent, both in am- 
plitude and phase, with the Compton-Getting prediction. 

IceCube is sensitive to muons from cosmic rays with 
TeV energies, and the data collected with this detector has 
been used to provide the first look at the CR anisotropy 
in the southern sky. The large set of cosmic ray events 
from IceCube, together with the air shower data from the 
IceTop detector, provide us with important tools to study 
the anisotropy of cosmic rays in the TeV and PeV range. 

The most recent results on CR anisotropy obtained 
with the IceCube and IceTop detectors are summarized 
in this paper. 

2. The IceCube and IceTop detectors 

IceCube is a km 3 neutrino telescope designed to search 
for astrophysical sources of high energy neutrinos. The ba- 
sic building block of IceCube is the Digital Optical Module 
(DOM), a glass pressure sphere that contains a 10" Hama- 
matsu PMT [14] , together with electronic boards for signal 
digitization [15], HV supply, and calibration LEDs. Be- 
tween 2004 and 2010, 5160 DOMs were deployed in the 
South Pole ice at depths between 1450 and 2450 m to de- 
tect the Cherenkov radiation emitted by charged particles 
produced in the interaction of neutrinos with nucleons. 

These DOMs are attached to 86 vertical strings that 
provide mechanical support, electrical power, and a data 
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connection to the surface. The vertical spacing between 
most consecutive DOMs in each string is about 17 m, while 
the horizontal spacing between most neighboring strings is 
approximately 125 m. 

A dedicated cosmic ray air shower array called IceTop 
is located on the ice surface above IceCube. The array 
consists of 81 stations, with two light-tight ice Cherenkov 
tanks per station. Each tank is 1.8 m in diameter, 1.3 m 
in height, and is instrumented with two DOMs that are 
operated at different PMT gains to increase the dynamic 
range of the detector. 

IceCube and IceTop were operated in partial configu- 
rations from the beginning of their construction until the 
completion of both detectors in December of 2010. Ice- 
Cube was operated in a 59-strings (IC59) configuration 
between May 2009 and May 2010, with IceTop operating 
with 59 stations (IT59) during the same time period. 

The IC59 dataset consists of those events where at 
least 8 IceCube DOMs detected photon hits within a 5 
lis window. The average rate for this trigger condition 
was 1.7 kHz in IC59. During a total live time of 334.5 
days, 3.4 x 10 10 events were collected, almost all of them 
produced by down-going muons from cosmic rays. For 
this analysis, a fast muon track reconstruction was per- 
formed online at the South Pole. The result of the fit, to- 
gether with the number of triggered DOMs and the time 
of the event are stored and transferred over a satellite 
link using a compressed data format. The median en- 
ergy of primary cosmic rays in this dataset is 20 TeV and 
was determined through Monte Carlo simulations assum- 
ing a mixed CR composition dictated by the polygonato 
model [16]. In this model, the energy spectrum for each 
chemical element is given by a broken power-law with a 
smooth transition, where the location of the spectral break 
is rigidity-dependent. Due to this dependence, heavier el- 
ements dominate the all-particle CR spectrum at energies 
above a few PeV. 

The median angular resolution of the muon track re- 
construction is 3°. Due to the degradation of the resolu- 
tion with increasing zenith angle, only events with < 65° 
were used in the analysis, which reduced the final dataset 
size to 3.2 xlO 10 events. 

In IceTop, the high-gain DOMs in the two tanks that 
form a station are run in local coincidence mode and the 
readout is enabled if they record hits within ±1 /is of each 
other. The IceTop trigger condition is satisfied if at least 6 
DOMs recorded locally-coincident hits within a time win- 
dow of 5 /is, which implies that at least 2 stations have 
participated in the event. 

The anisotropy analysis used events in which at least 3 
IceTop stations had triggered. Due to bandwidth limita- 
tions, events triggering less than 8 stations were prescaled 
by a factor of 8 while events with at least 8 stations were 
not prescaled. The event directional reconstruction was 
performed doing a x 2 fit to the trigger times of each station 
using a planar approximation for the shape of the shower 
front. Simulations show that the median resolution of this 



reconstruction algorithm is 2°. 

Preliminary results from Monte Carlo studies using a 
mixed composition of H, He, and Fe from the polygonato 
model indicate that the median primary CR energy of the 
IceTop dataset is 640 TeV, with 68% of the event between 
200 TeV and 2400 TeV. Only events with < 60° were 
selected for the analysis, with 1.2 x 10 8 events passing the 
cut. 

3. Analysis and results 

An anisotropy in the arrival direction of TeV cosmic 
rays was observed for the first time in IceCube using data 
from the 22-string configuration (IC22) that operated be- 
tween June 2007 and March 2008 and was reported in Ref. 
[T7] . In this analysis, the exposure-corrected right ascen- 
sion distribution of cosmic ray events was fitted with a 
harmonic function of the form Ai cos(i(a — fa)) + £>, 
where Ai and fa are the amplitude and phase of the i th 
term in the sum, a is the right ascension, and B is a con- 
stant. This sum was performed over the first two terms 
in harmonic space (n = 1,2) since they provide with an 
adequate description of the shape of the anisotropy. The 
fit parameters obtained in this analysis are A\ = (6.4 ± 
0.2 (sta) ± 0.8 (sys) ) x 10- 4 , fa = 66.4° ± 2.6° sta) ± 3.8° sys) , 
A 2 = (2.1 ± 0.3 (sta) ± 0.5 (sys) ) x 10- 4 , fa = -65.6° ± 

4 -°(sta) ± 7 ' 5 (sys) with X 2 / dof = 22/19, and show a good 
agreement with the phase and amplitude of the anisotropy 
observed in the northern sky. 

A later analysis [18] using IC59 data revealed that be- 
sides the large-scale structure (i.e. dipole and quadrupole 
modes) observed in the IC22 analysis there are also sta- 
tistically significant structures with typical sizes between 
10° and 20°. In this analysis, the search for anisotropy 
is conducted by searching for deviations of the sky map 
of reconstructed cosmic ray arrival directions in equatorial 
coordinates from a reference isotropic sky map obtained 
from data using the time-scrambling method described in 
Ref. [19] . The time scrambling period used in the analysis 
is 24 hours, which makes it sensitive to all angular scales 
in the celestial sphere. During the time scrambling proce- 
dure, events were resampled 20 times to reduce statistical 
fluctuations in the reference sky map. 

The sky maps were constructed using the HEALPbQ 
library [20] that provides an equal area pixelization of the 
sphere. The chosen HEALPix resolution divides the sphere 
into 49152 pixels, with an average distance between pixel 
centers of approximately 1°. Using the reference and data 
maps, a relative intensity map can be calculated using the 
expression 8U = (N { - (N)i)/(N)i, where iV* and (N)i 
are respectively the number of observed events and the 
number of reference events for the isotropic expectation in 
the i th pixel obtained with the time scrambling technique. 
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Figure 1: Angular power spectrum of the IC59 relative intensity 
skymap. The power spectrum before (blue) and after (red) the sub- 
traction of the dominant dipole (i = 1) and quadrupole (£ = 2) 
terms are shown for a time scrambling period of 24 hours. The me- 
dian value for the isotropic expectation is shown as a dashed black 
line, while the la and 2a bands are shown in gray. The power spec- 
trum shows significant departures from isotropy between £ ~ 4 and 
£ ~ 12 even after the subtraction, of the low-order terms. 



The angular power spectrum of the relative intensity 
map can be used to estimate the strength of the anisotropy 
at different angular scales in our data. The IC59 power 



spectrum is shown in Fig. 
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and was obtained using the 
that corrects for systematic 



effects introduced by the limited sky coverage of our data 
[2T] [22] . It can be seen that besides the already mentioned 
dipole (£ = 1) and quadrupole (£ = 2) there is a significant 
departure from isotropy at higher multipole moments be- 
tween £ ~ 6 and £ ~ 12, which corresponds to structures 
that have angular sizes roughly between 15° and 30° in the 
sky. In order to reveal this smaller structure, the dipole 
{£ = 1) and quadrupole (£ = 2) terms of the spherical har- 
monics functions were fit and subtracted from the IC59 
relative intensity map. The residual maps were smoothed 
to search for the small scale anisotropy. The smoothing 
procedure sums all events in a pixel to the events from 
pixels inside a certain angular distance. This produces a 
sky map of correlated pixels with an improved sensitiv- 
ity to structures with angular sizes similar to the smooth- 
ing radius. A scan was then performed over smoothing 
radii between 3° and 30° in 1° steps to find the optimal 
angular scale for the small scale structure. The optimal 
scale corresponds to the one that maximizes the statisti- 
cal significance of the observation calculated according to 
Ref. [23] , after taking into account trials due to the search 
over many pixels and smoothing radii. For the IC59 data, 
eight regions were identified where the absolute value of 
the statistical significance was higher than 5a before ac- 
counting for trials. These regions can be seen in Fig. [2] 
for smoothing radii of 15°. The most significant excess is 
region 1, with an optimal scale of 22° at which it reaches 
a significance of 5.3a after trials and has an amplitude of 
~ 10 -4 . A full list of the statistical significance for all 
regions can be found in Ref. [T8] , 
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Figure 2: Relative intensity sky map of the residual small scale struc- 
ture after the subtraction of a dipole and quadrupole terms for a 
smoothing radius of 15°. The labels correspond to the locations of 
all regions that showed an statistical significance larger than 5cr. See 
Ref. [18] for a full list of coordinates. 



A study of the energy dependence of the anisotropy 
was performed using IC59 data and was reported in Ref. 
[24] . A cut was implemented to create two distinct sub- 
samples with different median energies: 20 TeV, and 400 
TeV. The cut variables used in the selection were the re- 
constructed zenith angle of the event, and the number of 
triggered DOMs (both increase as a function of primary 
CR energy). After cuts, the 20 TeV dataset contained 
17.9 x 10 9 events, while the 400 TeV dataset consists of 
0.5 x 10 9 events (with 68% of the events between 100 TeV 
and 1300 TeV). The anisotropy results were obtained in 
two ways: through a harmonic fit to the right ascension 
distribution of events as in the case of IC22, and through 
the search for the optimal angular scale after a reference 
level estimation performed with the time scrambling tech- 
nique as in Ref. [18] . Both methods consistently observed 
the presence of the already known dipole and quadrupole 
structure in the 20 TeV dataset, while at 400 TeV the 
anisotropy pattern changes both in phase and amplitude. 
Only one structure in the 400 TeV sky map has a post- 
trial significance larger than 5a: a 6.3a deficit located at 
(a = 73.1°, S = -25.3°) with an optimal smoothing of 21° 
and an average amplitude of approximately 7 x 10 -4 in 
relative intensity. The deficit is also visible as region 6 
in the small scale map shown in Fig. [2] before energy cuts 
are applied. This observation represents the first detection 
of anisotropy in the southern sky at these energies. The 
relative intensity sky maps for the 20 TeV and 400 TeV 
energy bands are shown in Fig. [3] 

A preliminary analysis of the IceTop IT59 dataset re- 
veals a deficit located in the same region as the one ob- 
served at 400 TeV with IceCube. For a smoothing angle of 
20° the pre-trial significance is 6.2a and the amplitude is 
about 2 x 10 -3 , larger than the one observed in IceCube. 
A possible cause of this discrepancy is the difference in 
energy range associated with the two data sets. It is also 
possible that the CR chemical composition may be con- 
tributing to this mismatch. This is due to the fact that 
while IceTop is sensitive to all components of the CR air 
shower, IceCube is only capable of detecting the muon 
component, and this could create a detection bias towards 
a particular composition. Further studies of the energy 
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Figure 3: Equatorial sky maps for relative intensity in two different 
energy bands using the IC59 dataset: 20 TeV (above), and 400 TeV 
(below) for a smoothing radius of 20°. 



and composition dependence of the anisotropy are needed 
to perform a direct comparison between both results. 

A comparison of the right ascension projection of the 
relative intensities observed in IceCube for both the 20 
TeV and 400 TeV energy bands, and IceTop is shown in 
Fig. [5] Only events in the —75° < 5 < —30° declination 
range were used in this plot. 
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Figure 5: Relative intensity as a function of right ascension for the 
IC59 20 TeV (brown), and IC59 400 TeV (blue) compared to a pre- 
liminary result for the IceTop IT59 (orange) datasets. For clarity, 
only statistical error bars are shown. 



At energies of about 400 TeV, IceCube observed a 
strong deficit with a relative intensity of about 10 -3 and 
a size of approximately 20°. A preliminary analysis of Ice- 
Top data shows a deficit in the same region, but with an 
amplitude that doubles the one observed in IceCube. 

Future studies will expand the energy reach of the 
anisotropy analysis, and provide a handle on the evolu- 
tion of the anisotropy as a function of energy and angular 
scale. 
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Figure 4: Relative intensity sky map in equatorial coordinates for the 
IceTop IT59 dataset with a smoothing radius of 20° (preliminary). 



4. Summary 

Data taken between 2007 and 2010 with the IceCube 
neutrino telescope and the IceTop air shower array has 
been used to probe the anisotropy of TeV and PeV cosmic 
rays down to amplitudes of 10 -4 . The anisotropy at an 
energy of 20 TeV is consistent with that observed by other 
experiments in the northern hemisphere, and is dominated 
by a large scale component (dipole and quadrupole) with 
a strength of ~ 10 -3 . A subdominant, but statistically 
significant, structure at 20 TeV is characterized by small 
excess and deficit regions with angular sizes between 10° 
and 25° and strengths of the order of 10 -4 . 
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